Theory of the topological Anderson insulator 
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We present an effective medium theory that explains the disorder-induced transition into a phase 
of quantized conductance, discovered in computer simulations of HgTe quantum wells. It is the 
combination of a random potential and quadratic corrections cx p^a^ to the Dirac Hamiltonian that 
can drive an ordinary band insulator into a topological insulator (having an inverted band gap). We 
calculate the location of the phase boundary at weak disorder and show that it corresponds to the 
crossing of a band edge rather than a mobility edge. Our mechanism for the formation of a topo- 
logical Anderson insulator is generic, and would apply as well to three-dimensional semiconductors 
with strong spin-orbit coupling. 
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Topological insulators continue to surprise with unex- 
pected physical phenomena [1] . A recent surprise was the 
discovery of the topological Anderson insulator (TAI) by 
Li, Chu, Jain, and Shen 0. In computer simulations of 
a HgTe quantum well, these authors discovered in the 
phase diagram a transition from an ordinary insulating 
state (exponentially small conductance) to a state with 
a quantized conductance of Go = jh. The name TAI 
refers to the latter state. The findings of Ref. |2] were 
confirmed by independent simulations [5]. 

The phenomenology of the TAI is similar to that of 
the quantum spin Hall (QSH) effect, which is well under- 
stood lU [5] and observed experimentally in HgTe quan- 
tum wells [HI El [5] • The QSH effect is a band structure 
effect: It requires a quantum well with an inverted band 
gap, modeled by an effective Dirac Hamiltonian with a 
negative (socalled "topological") mass. The matching of 
this negative mass inside the system to the usual posi- 
tive mass outside leaves edge states in the gap. The edge 
states are "helical" , in the sense that the direction of 
propagation is tied to the electron spin. Opposite edges 
each contribute e^/ft. to the conductance. The conduc- 
tance remains quantized in the presence of (weak) dis- 
order, because time reversal symmetry forbids scattering 
between counter-propagating edge states (of opposite he- 

licity) UlSllSllZllH]. 

The crucial difference between the TAI and QSH 
phases is that the QSH phase extends down to zero dis- 
order, while the TAI phase has a boundary at a minimal 
disorder strength. Put differently, the helical edge states 
in the QSH phase exist in spite of disorder, while in the 
TAI phase they exist because of disorder. Note that the 
familiar quantum Hall effect is like the QSH effect in this 
respect: The edge states in the quantum Hall effect ex- 
ist already without disorder (although, unlike the QSH 
effect, they only form in a strong magnetic field). 

The computer simulations of Refs. [H [3] confront us, 
therefore, with a phenomenology without precedent: By 
what mechanism can disorder produce edge states with a 
quantized conductance? That is the question we answer 
in this paper. 



We start from the low-energy effective Hamiltonian of 
a HgTe quantum well, which has the form [5] 

H = a{p^a^-pyay) + {m + l3p^)a^ + [jp^ + U{r)]aQ. (1) 

This is a two-dimensional Dirac Hamiltonian (with mo- 
mentum operator p — —ihV, Pauli matrices a^, cTy, cr^, 
and a 2 X 2 unit matrix ao), acting on a pair of spin-orbit 
coupled degrees of freedom from conduction and valence 
bands. The complex conjugate H* acts on the opposite 
spin. We assume time reversal symmetry (no magnetic 
field or magnetic impurities) and neglect any coupling 
between the two spin blocks H and H* [S]. The scalar 
potential U accounts for the disorder. The parameters 
a,(3,j,m depend on the thickness and composition of 
the quantum well [7]. For the specific calculations that 
follow, we will use the same parameter values as in Ref. 
[2J, representative of a non-inverted HgTe/CdTe quan- 
tum well [TD]. 

The terms quadratic in momentum in Eq. ([l]) are not 
present in the Dirac Hamiltonian familiar from relativis- 
tic quantum mechanics, but they play an important role 
here. In particular, it is the relative sign of /3 and m that 
determines whether the clean quantum well {U = 0) is 
inverted (/3to < 0) or not-inverted (/3to > 0). We take 
/3 > 0, so the inverted quantum well has a negative topo- 
logical mass TO < 0. The inverted quantum well is a 
topological insulator (for Fermi energies Ep inside the 
gap), while the non- inverted quantum well is an ordi- 
nary band insulator. The phase transition between these 
two types of insulators therefore occurs at to = in a 
clean quantum well. 

We will now show that disorder can push the phase 
transition to positive values of to, which is the hallmark of 
a TAI. Qualitatively, the mechanism is as follows. Elastic 
scattering by a disorder potential causes states of definite 
momentum to decay exponentially as a function of space 
and time. The quadratic term Pp^ — —h'^pV^ in H, 
acting on the decaying state oc e~^^^, adds a negative 
correction dm to the topological mass. The renormalized 
mass TO — m + Sm can therefore have the opposite sign as 
the bare mass to. Topological mass renormalization by 



disorder, and the resulting change in the phase diagram, 
has previously been studied without the terms quadratic 
in momentum . The sign of m and m then remains 
the same and the TAI phase cannot appear. 

We extract the renormalized topological mass to, as 
well as the renormalized chemical potential fl, from the 
self-energy S of the disorder- averaged effective medium. 
To make contact with the computer simulations [2, ^3^ , we 
discretize iJ on a square lattice (lattice constant a) and 
take a random on-site disorder potential U, uniformly 
distributed in the interval {—Uo/2,Uo/2). We denote by 
-ffo(fe) the lattice Hamiltonian of the clean quantum well 
in momentum representation [51 [7] . 

The self-energy, defined by 
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with (• • • ) the disorder average, is a 2 x 2 matrix which 
we decompose into Pauli matrices: S — EqO'o + ^xCx + 
EyCTy -|- SzCr^. The renormalized topological mass and 
chemical potential are then given by 

TO^TO+limReS^, u = — lim Re En. (3) 

The phase boundary of the topological insulator is at 
fh — 0, while the Fermi level enters the (negative) band 
gap when \j2\ = —fh. 

In the selfconsistent Born approximation, E is given 
by the integral equation [T5] 

E = j^C/2(^/27r)2 / dk[EF + i(i^ -Ho{k)-Y.]-'^. (4) 

(The integral is over the first Brillouin zone.) The self- 
energy is independent of momentum and diagonal (so 
there is no renormalization of the parameters a,/?, 7). 
By calculating to and /I as a function of Ep and Uq we 
obtain the two curves A and B in Fig. [ij 

We have also derived an approximate solution in closed 
form |13j . 
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showing that the correction 6m = fh — m to the topo- 
logical mass by disorder is negative — provided /? > 7. 
For /? < 7 the clean HgTe quantum well would be a 
semimetal, lacking a gap in the entire Brillouin zone. 
Neither the TAI phase nor the QSH phase would then 
appear. In HgTe the parameter f3 is indeed larger than 
7, but not by much [TD]. Eq. ([5| implies that the lower 
branch of curve B (defined by /x = to < 0) is then fixed 
at Ep w TO > independent of Uq. This explains the 
puzzling absence of the TAI phase in the valence band 
{Ep < 0), observed in the computer simulations [Hl^- 
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FIG. 1: Computer simulation of a HgTe quantum well |10) . 
showing the average conductance (G) as a function of disorder 
strength Uo (logarithmic scale) and Fermi energy Ep, in a dis- 
ordered strip of width W = 100 a and length L = 400 a. The 
TAI phase is indicated. Curves A and B are the phase bound- 
aries resulting from the effective medium theory. Curve A sep- 
arates regions with positive and negative renormalized topo- 
logical mass fh, while curve B marks the crossing of the renor- 
malized chemical potential fl with the band edge {\p,\ = — fh). 
Both curves have been calculated without any adjustable pa- 
rameter. The phase boundary of the TAI at strong disorder 
is outside of the regime of validity of the effective medium 
theory. 



To quantitatively test the phase diagram resulting 
from the effective medium theory, we performed com- 
puter simulations similar to those reported in Refs. [Hl^- 
The conductance G is calculated from the lattice Hamil- 
tonian |SJ |7| in a strip geometry, using the method of 
recursive Green functions. The strip consists of a rectan- 
gular disordered region (width W, length L), connected 
to semi-infinite, heavily doped, clean leads [2]. Theory 
and simulation are compared in Figs, [l] and [2] 

Fig. [1] shows the phase diagram. The weak-disorder 
boundary of the TAI phase observed in the simulations 
is described quite well by the selfconsistent Born approx- 
imation (curve B) — without any adjustable parameter. 
Curve B limits the region where (A) the renormalized 
topological mass to is negative and (B) the renormalized 
chemical potential /2 lies inside the band gap: \p.\ < —fh. 
Condition (A) is needed for the existence of edge states 
with a quantized conductance. Condition (B) is not 
needed for an infinite system, because then Anderson lo- 
calization suppresses conductance via bulk states as well 
as coupling of edge states at opposite edges. In the rel- 
atively small systems accessible by computer simulation, 
the localization length for weak disorder remains larger 
than the system size (see later). Condition (B) is then 
needed to eliminate the bulk conductance and to decou- 
ple the edge states. 
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FIG. 2: Black curves, left axis: Average density of states p as 
a function of Ef for Uo = 100 meV, calculated by computer 
simulation (solid curve, for a disordered 100 a x 100 a square 
with periodic boundary conditions) or by effective medium 
theory (dashed curve). Red curve, right axis: Average con- 
ductance (G), calculated by computer simulation in a disor- 
dered strip (Uo = 100 meV, W = WOa, L = 400a). The TAI 
phase of quantized conductance lines up with the band gap. 
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FIG. 3: Red dashed curve: Average conductance (G) as a 
function of disorder strength {Ep = 25meV, W = L — 100 a). 
Black solid curve: Localization length ^, showing the peak 
at the strong-disorder edge of the conductance plateau — 
characteristic of a localization transition. The scaling with 
system size W of the width 5Uo of the peak is shown in the 
inset (double- logarithmic plot). 



Fig. [2] shows the average density of states p at the 
Fermi level. The agreement between the selfconsistent 
Born approximation (dashed black curve) and the com- 
puter simulation (solid black) is quite good, in particu- 
lar considering the fact that this plot is for a disorder 
strength which is an order of magnitude larger than the 
band gap. The range of Fermi energies over which the 
gap extends lines up nicely with the conductance plateau, 
shown in the same figure (red curve) . 

The strong-disorder phase boundary of the TAI cannot 
be described by effective medium theory, but it should be 
similar to the QSH phase boundary. In the QSH effect 
the strong-disorder transition is in the universality class 
of the quantum Hall effect JJj — in the absence of cou- 
pling between the spin blocks [121 [T7] . To ascertain the 
nature of the strong-disorder transition out of the TAI, 
we have calculated the critical exponent v governing the 
scaling of the localization length ^. For that purpose 
we roll up the strip into a cylinder, thereby eliminating 
the edge states [3]. We determine the localization length 
^ = —2 lim^^oo -Zj(ln G/Gq)"^ by increasing the length 
L of the cylinder at fixed circumference W. 

In Fig. [3] we show ^ as a function of disorder strength 
Uq at Ep = 25 meV, W = 100 a. As mentioned above, 
^ becomes much larger than W upon crossing the weak- 
disorder boundary of the TAI, so there is no localization 
there |18j . At the strong-disorder boundary, however, the 
dependence of ^ on Uq shows the characteristic peak of 
a localization transition [12]. In the inset we plot the 
scaling with W of the width SUq at half maximum of the 
peak. This yields the critical exponent via SUq oc W^^I'^ . 
We find V = 2.66 ± 0.15, consistent with the value v = 
2.59 expected for a phase transition in the quantum Hall 



effect universality class [20^. 

In conclusion, we have identified the mechanism for 
the appearance of a disorder-induced phase of quantized 
conductance in computer simulations of a HgTe quan- 
tum well "5^ . The combination of a random potential 
and quadratic momentum terms in the Dirac Hamilto- 
nian can change the sign of the topological mass, thereby 
transforming a non-inverted quantum well (without edge 
states in the band gap) into an inverted quantum well 
(with edge states). The weak-disorder boundary in the 
phase diagram of the TAI has been calculated by effective 
medium theory, in good agreement with the simulations 
(curve B in Fig. [T]). 

Contrary to what the name "topological Anderson in- 
sulator" might suggest, we have found that the hallmark 
of the TAI in the simulations, the weak-disorder tran- 
sition into a phase of quantized conductance, is not an 
Anderson transition at all. Instead, the weak-disorder 
boundary B marks the crossing of a band edge rather 
than a mobility edge. A mobility edge (similar to the 
QSH effect [15l [HI [17]) is crossed at strong disorder, as 
evidenced by the localization length scaling (Fig.[3|. 

Our findings can be summed up in one sentence: "A 
topological insulator wants to be topological". The 
mechanism for the conversion of an ordinary insula- 
tor into a topological insulator that we have discovered 
is generically applicable to narrow-band semiconductors 
with strong spin-orbit coupling (since these are described 
by a Dirac equation, which generically has quadratic mo- 
mentum terms [21j). There is no restriction to dimen- 
sionality. We expect, therefore, a significant extension of 
the class of known topological insulators [T] to disordered 
materials without intrinsic band inversion. 
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